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Towards multigrid methods for propagators of staggered fermions with 
improved averaging and interpolation operators 

T. Kalkreuter^* 

^11. Institut fiir Theoretische Physik, Universitat Hamburg, 
Luruper Chaussee 149, D-22 761 Hamburg, Germany 

A Dirac choice for the averaging kernel C is implemented numerically. This improved kernel will be needed in 
gauge covariant multigrid computations for propagators of staggered fermions. Results for C and the variational 
coarse grid operator will be given in 2-d SU{2) gauge fields. C-\ — |- is advocated for future algorithm development. 



Big efforts have been undertaken to find effi- 
cient multigrid (MG) methods for the computa- 
tion of propagators in gauge fields; see [1] for a 
list of references. Up to now no practical method 
has been found which is competitive for fermionic 
propagators in A-d non-abelian gauge fields. How- 
ever, it was proven that (geometric) MG works in 
principle without critical slowing down (CSD) in 
arbitrarily disordered systems, both for bosons 
[2,3] and for staggered fermions [4]. Unfortu- 
nately, this "ideal" algorithm is not practical be- 
cause the interpolation kernels have exponential 
tails whose truncation brings CSD back [3]. Nev- 
ertheless, the proof that CSD can be eliminated 
by MG raises hope that there exist also practical 
- though sophisticated - competitive methods. 

We consider the squared Dirac equation, 

i-p' + m')x = f , (1) 

where / may be a pseudofermion field, for in- 
stance. We will focus on staggered fermions and 
a gauge covariant "Dirac choice" [5] (no gauge 
fixing required) for the averaging kernel C . The 
C-kernel should map smooth functions on a fine 
grid onto functions which are smooth on a coarser 
scale. In [6,7] it was proposed to adopt an alge- 
braic definition of smoothness in disordered sys- 
tems. 

Given a differential operator D, according to 
[6,7] a function x is smooth on length scale a when 



11^X11 < 11x11 (in units a = 1) 



(2) 
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Figure 1. Norm of the ground-state of —p on 
the even sublattice of a 24^* lattice in a quenched 
SU(2) gauge field at /3 = 2.5. (Two-dimensional 
cut through the site with maximal norm.) 

Smooth functions m the sense of Eq. (2) are lo- 
calized m disordered systems. Baker [8] found in 
quenched 2-d SU (2) lattice gauge theory that the 
eigenfunctions corresponding to the 2 or 3 lowest 
eigenvalues of the negative gauge covariant Lapla- 
cian are localized for [3 = 0(1). A fermionic ex- 
ample in A-d SU(2) is given in Fig. 1. 

The localization property of the lowest mode of 
—p"^ could be taken into account in a geometric 
MG algorithm by shifting the site with maximal 
norm such that it coincides with a block center. 
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Let us now turn to the improved averaging ker- 
nel C . In Ref. [5] two qualitatively different pro- 
posals were made for the choice of C for stag- 
gered fermions in non-abelian gauge fields. The 
first one, the "Laplace choice" (LC), was used 
successfully in the "ideal" algorithm [4], but it 
has proven to be poor in practical variational MG 
methods [1,3]. Therefore we will test the second 
proposal which is called the "Dirac choice" (DC) 
of C. 

In the DC, as in the LC, C is defined by a 
ground-state projection method. This means that 
the adjoint C* of C fulfills a gauge covariant 
eigenvalue equation,^ 

{-lpl,D,.C*){z, x) = Mx) C*{z, x) , (3) 

where — ^ ^ is & block-local approximation 
of —ip'^ specified below, and \o{x) is its lowest 
(gauge invariant) eigenvalue, z denotes a site in a 
fine lattice, and a; is a coarse grid site ("a block"). 
C(x,z) is an element of the linear span of the 
gauge group. 

The solution of Eq. (3) is made unique by im- 
posing a normalization and a covariance condi- 
tion. The normalization condition reads^ 

l^J2\\C(x,z)\\^ = l (4) 

for all X, where Li, is the blocking factor (= 3), 
summation is over all fine grid sites within a block 
X, and the norm used is the trace norm. The 
covariance condition reads (for gauge group U(l) 
or 5(7(2)) 

Cix,x) = r{x)l (5) 

with X being the block center, and r(x) > 0. 

We use a blocking procedure for staggered 
fermions which is consistent with the lattice sym- 
metries of free fermions [5]. This forces us to 
choose Li, = 3. Even Li, are not allowed. 

Fig. 2 illustrates our choice of blocks in case 
of the DC. The different fermionic degrees of 

^ C depends on the gauge field although this is not indi- 
cated explicitly. 

^Because of the nontrivial overlap of blocks in the DC, we 
cannot retain the previous condition (7(7* = H. However, 
in pure gauges nothing changes. 



freedom are called "pseudofiavor" [5]. Different 
pseudofiavors are distinguished by different sym- 
bols in Fig. 2. Block centers x are encircled. They 
are spaced by which means that the coarser 
grid has times the number of sites of the 

finer grid. The boundaries of seven blocks x are 
marked. 
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Figure 2. Dirac-choice-blocking of a two-dimen- 
sional staggered lattice. 



The operator ip]^ d x (3) equals p'^ with the 
following boundary conditions (BC). For the sites 
which have the same pseudofiavor as the central 
site X we impose Neumann BC on the boundary 
of X, for all other boundary sites we use Dirichlet 
BC. Details can be found in [5,3]. 

The DC (3) has the following properties. In 
the limiting case of a pure gauge, we regain the 
LC with nonoverlapping blocks. The C-kernel 
for free fermions is piecewise constant on sublat- 
tices of given pseudofiavor.^ In nontrivial gauge 
fields - contrarily to the LC - the averaging kernel 
C{x, z) of the DC is nonvanishing also when x and 

■^In pure gauges the variational MG algorithm with these 
kernels is successful in eliminating CSD [3]. 
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z carry different pseudoffavor (but have the same 
parity, i. e. are both even or both odd). Therefore 
the blocks in Fig. 2 overlap in a nontrivial way. 
This overlap takes the field strength term (T^^i^^y 
in Jj)'^ into account. (T^^i^^y had been disregarded 
in the LC. 

The DC for C has been implemented in C++ 
[9,10]. The author advocates this language for 
future algorithm development. It requires some 
effort to define suitable C++ classes, but once 
this is done, all advantages of C++ [10, Sec. 1] 
are at your disposal. 

The eigenvalue equation (3) was solved by the 
efficient algorithm of [11]. All computer programs 
were checked for gauge covariance (i. e. without 
gauge fixing corresponding results were found for 
different points on gauge orbits). 

First tests have been performed in quenched 
%d SU(2) gauge fields. An encouraging result 
is that Xo(x) = for all blocks x and for arbi- 
trary [3. Global lowest eigenvalues are collected 
in Table 1. Examples for C, for the variational 
coarse grid operator —Cip'^C*, and for the non- 
diagonal(!) mass term CC* are in Table 2. Fur- 
ther results in c? = 2,4, including computations 
of propagators, will be available shortly [12]. 

I am indebted to G. Mack for many stimulating 
discussions. Financial support by Deutsche For- 
schungsgemeinschaft is gratefully acknowledged. 
I wish to thank DESY for providing resources on 
HP workstations . 

Table 1 



Lowest eigenvalues; cf. Table 1 of [4]. 





|A| 


-_^2 


-cp'c* 


CXD 


any 








40 


242 


0.004513 


0.206828 


1.0 


242 


0.000173 


0.418796 
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Table 2 

Examples of kernels in 2-rf 5'[/(2) at /3 = 40 (^ f« 
5.2). The bold entries vanish in pure gauges. 





0.997 





0.996 





0.993 





0.052 





0.203 





1.00 
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1.10 
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0.984 
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3||- 


Cp'C*{Q,x)\\ 
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3.94 
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0.913 








ncc*{(),x)\\ 




0.060 





0.130 
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0.069 





0.078 
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